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1. INTRODUCTION 



The physics of QCD exhibits different relevant excitations at distinct length (or momen- 
tum) scales. To give this notion a well-defined meaning, consider some characteristic length 
scale L c of the order of 1 fm that crudely separates short- from long-distance physics. At 
short space-time distances (r <C L c ) the relevant degrees of freedom are quarks and gluons, 
effectively unconfined due to asymptotic freedom, and their interactions are well described 
by perturbative QCD. The theory exhibits chiral symmetry and (approximate) scale sym- 
metry. At large distances (r S> L c ) on the other hand, we are in the regime of hadronic 
degrees of freedom and physical observable particles, whose non-perturbative interactions 
are known to be described well by chiral models. In between these two regimes, in the 
range r ~ L c , our current knowledge is essentially limited to the understanding that there 
must be a rather sudden dynamical establishment of long-range order, i.e. some kind of 
"phase transition" from the unconfined, chiral- and scale-invariant phase of partons to the 
hadronic phase with massive physical states and broken symmetries. 

The dynamics of this parton-hadron conversion and confinement mechanism has scarcely 
been studied yet, although QCD-inspired effective quark models that incorporate confine- 
ment phenomenology in some way have been exploited extensively to describe static hadron 
properties rather well M. This problem is particularly serious for attempts to describe the 
phenomenon of hadronization in high-energy QCD processes. The theoretical tools cur- 
rently available for studying QCD are inadequate to describe the transformation from par- 
tonic to hadronic degrees of freedom as a dynamic process: perturbative techniques j|] are 
limited to the short-distance regime where confinement is not apparent, whilst effective low- 
energy chiral models || and QCD sum rules ||], that incorporate confinement, lack partonic 
degrees of freedom. On the other hand, common descriptions of parton fragmentation ||] 
are usually based on ad hoc prescriptions to simulate hadron formation from parton decays. 
In principle, lattice QCD Q should be able to bridge the gap, but in practice dynamical 
calculations of parton-hadron conversion are not yet feasible. 

The purpose of the present paper is to give a detailed documenation of our progress 
towards a consistent, fully dynamical formulation of the non-static properties of confine- 
ment, chiral symmetry breaking and hadron formation, as recently proposed in Ref. j?]]. 
Aside from the aforementioned arguments, these issues are of great interest in the context 
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of the QCD phase transition in the early Universe when hadrons formed from unconfined 
quark-gluon matter, or in high-energy heavy-ion collisions, where one expects a very hot 
and dense deconfined quark-gluon plasma to be created. Here it is inevitable to employ a 
dynamical treatment of the transition from short-distance (perturbative) regime of partons 
to the long-distance (non-perturbative) domain of hadrons. We extend here previous work 
H and present a universal approach to the dynamic transition between partons and hadrons 
based on an effective QCD field theory description and relativistic kinetic theory. 

Our concept is the following: we start from a gauge-invariant Lagrangian formulation 
that embodies both fundamental partonic and composite hadronic degrees of freedom. It 
is explicitly dependent on the space-time scale L = \fr^ at which the physics is "probed" . 
The scale dependence is however not external, but the variation of the scale is governed by 
the dynamical evolution of the physical system under consideration. The field equations of 
motion can be cast into evolution equations for the real-time Green functions of the various 
particle species, and by following the space-time evolution we can trace the conversion from 
partonic to hadronic degrees of freedom in full 7-dimensional phase space, as it is driven 
by the dynamics. This effective field theory approach recovers QCD with its scale and 
chiral symmetry properties at short distances or high momentum transfers, but yields at 
low energies the formation of symmetry breaking gluon and quark condensates including 
excitations that represent the physical hadrons. 

It is important to explain in more detail the physical basis of our approach: we as- 
sume that the vacuum state in QCD can be visualized as a "color dielectric medium" ]|] 
chracterized by some collective color-singlet fields that correspond in the long-wavelength 
limit to the gluon and quark condensates, and incorporate phenomenologically the complex 
structure of the physical vacuum as order parameters. Specifically, the underlying hypoth- 



ysis [|1C], 11 1 is that the long-distance (non-perturbative) gluon self-interactions generate an 



effective scalar gluon condensate field x which is self-interacting through some potential V 
constructed [fLl] ] on the basis of the symmetry properties of the QCD Lagrangian. As a 
consequence of symmetry constraints, the scalar field x must in addition couple through 
the potential V to an effective pseudoscalar quark condensate field U, a feature which is 
also suggested by lattice QCD studies Q, indicating that the confining and chiral symme- 
try breaking "phase transitions" are in some way related and occur approximately at the 
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same scale. The long-range properties of the non-perturbative vacuum are then character- 
ized by gluon and quark condensates which are proportional to the non-vanishing vacuum 
expectation values (x) = Xo an d (U) = Uq in the long-distance limit. 

Our central idea is that the effective gluon condensate field x plays the driving role in the 
generation of confinement and chiral-symmetry-breaking mechanisms: the non-perturbative 
gluon self-interactions are assumed to modify the long-range properties of the vacuum in 
such a way that the propagation of the elementary gluon fields is altered with increas- 
ing space-time distances and eventually completely suppressed (confinement). As a direct 
consequence [fUJl , the self-energy of the elementary quark fields ijj will be modified accord- 
ingly through the quark-gluon coupling, so that it generates dynamically an effective quark 
mass term which becomes infinite at large space-time distances. The coupling between the 
perturbative regime with elementary fields A^, if) and the non-perturbative vacuum repre- 
sented by the condensate fields x, U is mediated by a single dimensionless "color-dielectric 
function" k(x) which vanishes in the short-distance limit, but approaches unity at 
large distances. There is no need to introduce an additional coupling to the field U, nor 
to consider an explicit chiral-symmetry-breaking quark mass term, because, as mentioned, 
the quarks aquire a dynamic mass via their coupling to the gluons in the presence of the 
field x- 

We thus obtain an effective QCD field theory which in the short distance limit ((x) = 
0, (U) = 0, k(0) = 1) is chiral invariant and incorporates free gluon and quark propagation 
(asymptotic freedom), whereas in the long distance limit ((%) = xo> (U) = Uo,k(xo) = 0) 
no gluon or quark propagation can occur (confinement). In between these two regimes, the 
effective theory interpolates and governs the dynamics of the conversion of short-distance 
fluctuations (partons) to non-perturbative bound states (hadrons) embedded in the physical 
vacuum. As a prototype case, we study in detail the parton-hadron conversion in e + e~ — > 
hadrons. We visualize the process e + e~ —* qq as producing a "hot spot" in which the long- 
range order represented by xo and Uo is disrupted locally by the appearance of a bubble 
of the naive perturbative vacuum. Within this bubble, a parton shower develops in the 
usual perturbative way, with the hot spot expanding and cooling in an irregular stochastic 
manner described by QCD transport equations. This perturbative description remains 
appropriate in any phase-space region of the shower where the local energy density is large 
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compared with the difference in energy density between the perturbative partonic and the 
non-perturbative hadronic vacua. When this condition is no longer satisfied, a bubble of 
hadronic vacuum may be formed with a probability determined by statistical-mechanical 
considerations. 

This paper is organized in two main parts. The first part, consisting of Sees. 2 and 
3, is intended as a comprehensive presentation of the general field-theoretical framework 
and necessary elements of quantum transport theory. In Sec. 2 we construct on the basis 
of the dual vacuum picture of coexisting perturbative and non-perturbative domains an 
effective theory that embodies the correct scale and symmetry properties of QCD and that 
has the desired features outlined above. We also discuss the relation to the phenomenology 
of the QCD phase transition, where the role of the critical temperature is analogous to the 
critical confinement length scale in our approach. Sec. 3 outlines the method of real-time 
Green functions that we use to derive from the field equations of motion the corresponding 
coupled equations for the particle distribution functions. We also indicate how macroscopic 
quantities related to observables can be extracted from the microscopic particle dynamics 
within the kinetic theory of (non-equilibrium) many-particle systems. The second part 
of the paper, Sees. 4-6, is devoted to the application of this effective QCD field theory 
to the dynamics of parton-hadron conversion for the prototype process of fragmenting jet 
systems initiated by e + e~ annihilation. In Sec. 4 we derive transport equations that, in 
the case of the partons, are generalized QCD evolution equations in full phase space, and 
similar equations for the excitations of the x an d U fields. In Sec. 5 we present results of 
simulating this real-time evolution of partons through the perturbative shower stage, via 
subsequent formation of color-singlet clusters, and finally hadronic cluster decay to give the 
final hadron yield. We investigate phenomenological implications for particle production 
and the Bose-Einstein effect, which we find to be a particularly sensitive probe to measure 
and test the confinement dynamics. Finally, Sec. 6 is reserved for a summary and a brief 
discussion of future perspectives of the approach, in particular its applicability to the QCD 
phase transition, and to high-density QCD. 
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2. FIELD THEORY FRAMEWORK 



As explained in the introduction, the vacuum state in QCD may be pictured as a color- 
dielectric medium characterized by long-range order parameters. Consider the vacuum 
expectation values (uev's) of the normal-ordered products of color-singlet and Lorentz scalar 
(pseudoscalar) functions x (U) [Uj 

(0| X (FD|0) ; <0|£/(^,^-)|0>, (1) 

where F£ v is the usual SU(3) field strength tensor, ip, tp the quark fields, and x(0), U(0) is 
set to be zero. For instance, x can De fabcF^u^uX^Xfi or {^v^v) 2 ■> or other combinations. 
Similarly, U may be composed of Trf^i^j) 2 ], Tr\{ij) i 'ij}jij}^if , \, etc.. These vev's are phys- 



ical quantities that characterize the structure of the QCD vacuum [16] and are related to 
the measurable gluon and quark condensates, respectively. 

Clearly, if we take the long wavelength limit L — > oo and simultanously let the coupling 
strength g s among the fields tend to zero, we have in general the non-comutativity of the 
double limits 

lim lim ( | X (FD |0) = + lim lim ( | X (^D I ) 

L— >oo g s — »U gs— >U L^oo 

lim lim ( | Utyi,^) I 0) =0 / lim lim ( 1 Utyi^A | 0) , (2) 

with (0|F^|0) = = (0|V>i|0). Eq. (Q) is a pure quantum phenomenon and a typical 
property of phase transitions. It implies that there is long-range order in the vacuum which 
can be characterized by the operator functions x an d U. 

In order to embody this concept into a field theory formulation, let us define the distance 
measure L for the space-time separation between two points r and r' (r** = (t,r)): 



L 



^{r -r>)^{r -r'Y , (3) 



and introduce a characteristic length scale L c that separates short distance (L <C L c ) and 
long range (L>L C ) physics in QCD. The scale L c can be associated with the confinement 
length of the order of a hadron radius, as we will specify more precisely later. 
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2.1 The short-distance regime L <C L c 

At small space-time distances L <C L c , because of asymptotic freedom, the properties 
of QCD are well described by a perturbative expansion in powers of the coupling g s of the 
generating functional for the connected Green functions, 

W[J,r],7]] = j VA»Vxl)V$ ds&F exp |t J d 4 r(£[^,V,^] + J^aK + ^iV + #»)} ■ 

(4) 

In the path integral, detj 7 denotes the Fadeev-Popov determinant and J, rj, rj are the 
generating currents for the gluon fields A 1 ^ and the quark fields if), tp (which are vectors in 
flavor space, ip = (V'uiV'rf; • • ■))> an( ^ the QCD Lagrangian is 



ll>3 + , (5) 



where F{f u = d^A v a — d u A^ + gsfabc^b^c 1S the gluon field-strength tensor. The subscripts 
a, b, c label the color components of the gluon fields, and g s denotes the color charge related 
to a s = <7g/(47r). The T a are the generators of the SU(3) color group, satisfying [T a ,Tj,] = 
ifabcTc with the structure constants f a bc- The indices i,j label the color components of 
the quark fields and m = diag(m M , ra^, . . .). Throughout, summation over the color indices 
a, b, c and i,j is understood. We recall that on setting the quark current masses m to zero, 
one has exact chiral symmetry. The gauge-fixing term is denoted by a general function 
£ a (^4) which, e.g., in covariant gauges is £ a (^4) — — l/(2a)(c^^4^) 2 with Lagrange multiplier 
1/a. However, we will later consider a different (ghost-free) gauge that is more convenient 
for our purposes. 

2.2 The long-distance regime L ^> L c 

The long-range physics of QCD at large space-time distances L S> L c , is known to be 
described well by an effective low-energy theory. Here we adopt the approach of Ref. [JTJ] 
and define the corresponding generating functional as: 



W 



■J x ,Ku,Kl] = J V X VUVrf expji J d A r(c{x,U,U^\ + J xX + U^K V + K^U 



(6) 

The field degrees of freedom are a scalar gluon condensate field \ and a pseudoscalar quark 
condensate field U = f n exp (i J2^=o ^j4>j/ fir) f° r the nonet of the meson fields <pj (f n = 93 
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MeV, Tr[\i\j] = 25-ij, UU^ = /^), with non- vanishing vev's in the long-distance limit, 



Xo 



Un 



<5J, 



-lnW[J x ,Ku,Kl] = i -^y L 
+ J-)lnW[J x ,K u ,Kl] -- 



(o| u + tf |o; 
<o|o> 



+ o 



(7) 

(8) 



and an effective action 



T[ X ,U,U^] = hxW[J x ,Ku,Kl] - J d 4 r{j xX + U^K V + K^u} 



0) 



+ 



Consequently the Lagrangian in (|6|) is given by 

1 



v( x ,u) , 



(10) 



with a potential V that has been constructed [12, 13] on the basis of constraints which arise 
from the scale and chiral symmetry properties of the excact QCD Lagrangian, namely, 



V(x,U) 



4 i 4 l 

4X0 + X ln Vpl/4 



X 



e^Xo 



1 

+ 4 



1 



2L\ 

Xo) 



Tr 



+ cTr 



rh q (U + U^) ( 



X_ 

.Xo 



1 2 ,2 f X 
z VXo 



(11) 



Here the parameter b is related to the conventional bag constant B by 

B 



b^. 



(12) 



Furthermore, c is a constant of mass dimension 3, m q = diag(m u , rrid, m s ) is the light 
quark mass matrix, and nig is an extra U(l)-breaking mass term for the ninth pseudoscalar 
meson <po (which we will disregard in the following). In the chiral limit, this potential has 
a minimum when (x) = Xo and equals the vacuum pressure B at (x) = 0. 

2.3 The intermediate regime L w L c 

Having established a field theory framework for the two regions L <C L c and L S> L c , the 
crucial issue is now the intermediate range. Clearly there must be a dynamical interpolation 
around L ~ L c from the short-range to the long-range description. We propose here the 
following approach. Let us first consider the long-range domain, i.e. the physical vacuum 
characterized by Xq, and introduce into the vacuum an excitation of small space-time extent. 
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For instance, imagine the creation of a qq pair with invariant mass Q ~ L _1 3> L^ 1 
by a time-like virtual photon from e + e~ annihilation. The insertion of such a localized 
excitation ("hot spot") modifies the vacuum and we assume that the corresponding change 
in the action integral S = J d^rClx, U, U>] in (B) can be evaluated sufficiently accurately to 
second order as 

5 2 C[ X , U, l/t] 



5S 



d l r 



) F^ a F b ^, 



SF^aSFf 



+ o 



5 2 C[x,U,U^} 



) ipitpj 8ij J 
(13) 



where k_l and /i£ refer to the appropriate feu's. Note that this change 5S in the action 
preserves local gauge invariance. We also remark that this ansatz implicitly assumes that 
the elementary gluon (F^ v ) and quark fields (t(),ip) couple directly only to the scalar field 
X, but not to the pseudoscalar field U. The dynamics of U is solely driven by its coupling 
to x through the potential V(x, U), eq. (|PI|). 

On the other hand, we know that the short-range properties at L -C L c of our qq 
excitation are not affected by the long-range correlations. Thus, here we can use (|) with 
perturbative methods, since the quanta are asymptotically free and (x) = = (U + U'). 
Thus we can combine @ and the effect of (13) by adding to C[A^, i/j, and jO[x,U,W) 
the following contribution that carries an explicit scale- (L-)dependence: 



C L [A",il,,il>,x] 



1 - kl(x) 



- kl(x)) UA) 

(14) 

where the third term in the integrand is necessary to maintain local gauge invariance. It 
remains to specify the form of the functions kl(x) an d Hl{x)- Since kl has to satisfy the 
boundary conditions |l5j 



«i(0) 



1 , 



0, 



(15) 



and is constrained to be a Lorentz-invariant color-singlet function of scale dimension zero, 
a minimal possibility is 

«fc(x) = 1 " (t^ttY ■ ( 16 ) 



Lq Xo 



It turns out that the particular form of kl(x) ls n °t crucial as long as the properties (15) are 



satisfied fll7 ], because parton-hadron conversion is quite rapid, as we see later, being related 
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to the weakly first-order nature of the QCD phase transition at finite temperature. The 



essence is that (16) enforces color charge confinement due to the fact that a color electric 
charge creates a displacement D a = K L E a , where E% = F® k , with energy \ J d 3 rDl/K L 
which becomes infinite at large r for non-zero total charge. 

Similarly, absolute confinement can be ensured also for quarks by coupling the quark 
fields to the x field through 

" l(x) = m Uzw - V = (w - (L X r • <17) 

where hq is a constant of mass dimension one that we will set equal to 1 GeV. This form 
reflects that the quark mass term Hl{x) in (0) I s induced by non-perturbative gluon 
interactions, rather than being an independent quantity, as is suggested by an explicit 
calculation [11] of the quark self-energy involving the gluon propagator in the presence of 



the collective field x- It h as been shown [ 10 1 that the dynamical mass hl{ X ) leads to an 



effective confinement potential with the masses of the quarks at small L approximately equal 
to the current masses, but at large L when (x) — > Xo it generates an infinite asymptotic 
quark mass, 

Hl(0) = , Hl(xo) = oo . (18) 



It is evident from (|iq)-(|l8|) that C[A^, if), tp, x] given by (14) vanishes in the short-distance 
limit (L — > 0, (x) ~* 0)) whereas in the long-distance limit it suppresses the propagation of 
colored gluon and quark fluctuations, and interpolates smoothly between the two extremes. 
The typical functional forms of kl(x) an d A*Z/(x) are illustrated in Fig. 1. 

2.4 The scale-dependent generating functional for the effective theory 

Let us now summarize and combine the three contributions of Sees. 2.1-2.3 into a 
single action integral, and write down the resulting generating functional as an effective 
description covering the full range < L < oo and depending implicitly on the scale L as 
defined by @: 

W L [J,r},r},J x ,K u ,K'} / ] = J VA^V^VxVUVU^ detF (19) 

x expji J d 4 r(£[^,V,?] + C L [A»,1)rf, X ] + C\x,U,rf] 



+ J^A* + iPv + fj^J xX + U^Ku + Kpty 



where C[A», ip, tp] is given by (§), C[x, U, U^} by (0), and C L [A^, tp, ip, x] b Y ©, so that the 
effective, L-dependent Lagrangian density £^ = £[^4/*, tp, tp] + C[x, U, W] + Cl[A^, tp, tp, x] 



in the generating functional (19) can be written as 



~(^x)(^x) + ~ 



(cL£/)(cWt) 



^7 + «i(xKa(^) 



where V(x,U) is the potential given by (11), and we have gone over to the limit of zero 



current quark masses. It is important to realize that the scale dependence of (19) and 
( f20l) arises solely through the L-dependent functions kl(x) an d Ml(x); given by eqs. (lilf ) 
and (|l^). The scale L is not to be misunderstood as an external parameter. Instead it is 
intrinsic variable of the formulation. As we will see later, the variation of L is governed 
by the dynamics of the fields itself, and it in turn determines the time evolution of the 
interacting fields. Therefore, when studying the dynamical evolution of some system under 
consideration, one must necessarily require this self-consistency for a meaningful solution. 

At this point let us state clearly the following important remarks: 

a) The effective field theory defined by ([l9|) and (|20|) represents a description of the 
duality of partonic and hadronic degrees of freedom: high-momentum, short-distance quark- 
gluon fluctuations (the perturbative excitations) are embedded in a collective field x (the 
non-perturbative vacuum), in which by definition the low- momentum, long-range fluctu- 
ations are absorbed. Confinement is thus associated with a dual stucture of the QCD 
vacuum. The formulation is gauge- and Lorentz-invariant, and is consistent with scale and 
chiral symmetry properties of QCD. It interpolates between the high-momentum (short- 
distance) QCD phase with unconfined gluon and quark degrees of freedom and chiral sym- 
metry ((x) = 0, (U) = 0,kl = 1,/^l = 0), to a low-energy (long-range) QCD phase with 
confinement and broken chiral symmetry ((%) = xo, (U) = Uq, kl = = oo), where xo 
and Uq are the long-range order parameters of the vacuum, directly related to the gluon 
condensate and the quark condensate, respectively (Sec. 2.5 below). 

b) By introducing additional fields x an d U to describe the long-range behaviour of the 
gluon and quark fields, we must obviously be careful not to double-count the degrees of 
freedom, since the full theory of QCD is contained in £[y4/*, if), ip], eq. (||), already. However, 
by our construction the sum Cl in (|2^) gives a consistent formulation that strictly avoids 
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double counting, because the introduction of the scale L and the behaviour of the In- 
dependent coupling functions k_l(x) an d Hl{x) truncate the dynamics of the elementary 
fields A^, to the short-distance, high-momentum regime (L <C L c ), whereas the effective 
description in terms of the collective fields Xi U covers the complementary long-range, low- 
energy domain (L 3> L c ). Accordingly, a quark or gluon is either considered a colored 
short-range fluctuation (parton) or it is part of a complex bound state (hadron), but not 
both. 

c) The presence of the non- linear coupling function hl{x)i which also enters (J,l{x) v * a d0~ 



means that the sum AC = C[x, U] + C[A^, tp, ip, x] in (20) is non-renormalizable. However, 
there is no need for explicit renormalization, because the composite fields x an d U are 
already interpreted as effective degrees of freedom with loop corrections implicitly included 
in AC, and it would be double counting to add them again. Moreover, as mentioned in 
item b) above, the low-energy domain of C[x, U] is by construction bounded from above 
by the onset of the high-energy regime described by C[A^, tjj, ip]. The characteristic scale 
L c that separates the two domains, therefore, provides an 'ultra-violet' cut-off for C[x, U], 
and at the same time an 'infra-red' cut-off for £[A^, ?/>, 

2.5 Analogies with QCD at finite temperature 

We close this Section with pointing out some immediate phenomenological implications: 
the particular form ( |TT| ) of the potential V as a function of x> as well as the functions kl and 
jiL that couple short- and long-range regimes, play a central role in dynamical processes 



where the scale L changes with time. The effect of (14) can be interpreted as a scale- 



independent modification 5V , which adds to the (L-independent) potential V, eq. (|l l|) , 

V(L) := V( X ,U) + 5V(L, X ) , (21) 

™ = W F ^ Fr + (W-1lx)» * ' (22) 

where we used eqs. ( |l6| ) and (p^7[). We emphasize again that L is not an input parameter, 
but rather is determined by the space-time dependent separation of the colored quanta. In 
Sec. 4, we will specify how to determine the variation of the variable L. 

In view of (p2[), one has SV oc 0(L 2 ), therefore it is suggestive that the variable L 
plays a similar role as the temperature T in finite-temperature QCD, where the correction 
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to the zero-temperature potential is 0(T 2 ) |18[] . This formal analogy will be indicative in 
the following. However, one must bear in mind that here we are aiming to describe the 
evolution of a general non-equilibrium system in real time and Minkowski space, as opposed 
to thermal evolution in Euclidean space. Nevertheless, we adopt the general concept of Ref. 
jnj, and from the analogy with this previous work we can qualitatively expect that the 
correction 5V will give a first-order "phase transition" from the parton to the hadron phase, 
when combined with V. 

As seen in Fig. 2, there are three characteristic scales, L x , L c and Lo, that mark the 
time evolution from the small- L to the large-L region as the scale-dependent potential V(L), 
eq. (|2"l|), changes: 

(i) L x is the characteristic length scale below which the vacuum with x / cannot exist. 
The potential V has a unique minimum at x = 0, i.e. we are in the perturbative vacuum 
of the pure parton phase. 

(ii) L c marks the point when the V develops two degenerate minima, one at x = an d 
the other at \ = Xc- The pressure in the parton phase is here equal to the pressure in the 
hadron phase, and the probability for partons to tunnel through the barrier becomes large. 



(iii) Lq defines when 5V = and V becomes equal to V in eq. (|llj), and has a single 
absolute minimum at (x) = Xo- The parton phase cannot exist any longer, and the parton- 
hadron conversion is completed. We are in the true (physical) vacuum characterized by the 
presence of a gluon and a quark condensate. 

Following we can relate the vev's (!§), xo and Uq, to the gluon condensate 

( 0|^i^FH0> = -bxt = Go (23) 
and the quark condensate 

<0|gg|0) = c (Z-) 3 U = Qo, (24) 

where b and c are defined in (|j~l|). These condensates can be regarded as local order pa- 
rameters associated with gluon and quark confinement, respectively, and chiral symmetry 



breakdown. Also, as discussed in [14], one can interpret small oscillations about the mini- 



mum of the potential V = V at (x) = Xo> (U + U'} = Uq, as physical hadronic states that 
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emerge after symmetry breaking. They include in principle: (a) glueballs and hybrids as 
quantum fluctuations in the gluon condensate xoi (b) pseudoscalar mesons as excitations 
of the quark condensate Uq, (c) the pseudoscalar flavor singlet meson fa, and (d) baryons 
as topological solitons. We will return below to the issue of hadron formation. 

3. EQUATIONS OF MOTION AND KINETIC EVOLUTION 

In this Section we outline how to obtain a fully dynamical description of the parton- 
hadron conversion in real time and complete phase space, starting from the defining gener- 
ating functional (|i~9| ) of our effective field theory. A comprehensive derivation can be found 
in Ref. to which we refer for details. The method is based on the Green function 



technique p0| , 21], here applied to derive transport equations for the field operators A 11 , ip, 
X, and U. The form of the transport equations results directly from the Dyson- Schwinger 
equations The self-energy operators that enter the connected part of the equations 

can then be evaluated in a perturbative expansion. This leads to corresponding equations 
of motion for the distribution functions of particles, namely gluons and quarks as colored 
fluctuations, and scalar and pseudoscalar hadronic excitations. The solution for the time 
development of these particle distribution functions in phase space will then allow us to cal- 
culate macroscopic observable quantities within the framework of relativistic kinetic theory. 



3.1 Equations of motion 

We start from the field equations of motion that follow from the variation of the gener- 
ating functional © with @: 

1 1>i = 



fy*T + 9sf abc A,, b Fr - (a M ln KL (x)) Ff ~ (l + ^) 9s ^ + Q(A) 



where = d/dx^, and we set U = (U + W)/2. In the second equation, the function 

dUA) 



C(A) := dAK L §^L\ (26) 
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results from the gauge-fixing constraint in eq. (||), e.g., in covariant gauges £„(.A) = 
(l/a)d u duA£. Note that in general there are additional equations of motion involving 
the ghost fields coupled to the gluon fields, however, since we will later choose a ghost free 
gauge, these decouple and are irrelevant here. 

It is evident from the above equations and the form of the potential V(x, U), eq. (|TT|) , 
that in the short-distance limit when (x) = and k>l{x) = 1> the system of equations 
decouples and reduces to the usual Yang-Mills equations. Similarly, in the long-wavelength 
limit, one has (x) — ¥ Xo, (U) — > Uq, and kl{x) ~* 0> so that the dynamics in this case is 
completely described by the equations for the effective fields x an d U. A very important 
point is that the U field does not couple directly to the quark or gluon fields. By construction 
[14 1, the dynamics of the quark condensate field U is solely driven by the gluon condensate 



field x- As a consequence, the equation for U is readily solved, once the solution for x is 
known. It is important to realize that the interplay between the x field and the quark and 
gluon fields, if) and A, is the crucial element of this approach. From a phenomenological 
point of view, this implies that the transformation of parton to hadron degrees of freedom 
proceeds first by formation of scalar color-singlet states which subsequently decay into 
pseudoscalar excitations (Sec. 4). 

3.2 Real-time Green functions and microscopic kinetics 

The central role in the following is played by the real-time Green functions in the 'closed- 
time-path formalism' (for an extensive review, see pl| ). This formalism is the appropriate 
tool to describe general non-equilibrium systems, and its particular strength lies in the 
possibility of studying the time evolution of phenomena where initial and final states cor- 
respond to different vacua, as we are addressing here. In |l!| it is shown how one obtains 
a dynamical formulation which systematically incorporates quantum correlations and de- 
scribes naturally the transition from the perturbative QCD regime to the non-perturbative 
QCD vacuum. The real-time Green functions are defined as the two-point functions that 
measure the time-ordered correlations between the fields at space-time points x and y (as 
before we suppress the spinor indices for the fermion operators): 

iSufay) := (Tijjiix^M) 



'ah 



xD%{x,v) := (TA^x)AUy)) 
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iA(x,y) := {Tx{x)x{y)) ~ ( xO) } ( x(v) ) > 

i A(x, y) := (T U(x) U^y) ) - ( U(x) )(U\y)), (27) 



where (. . .) denotes the vacuum expectation value, or, in a medium, the appropriate en- 
semble average, and T A(x)B(y) = 6(xo,yo)A(x)B(y) db 0(yo, xo)B(y)A(x) stands for the 
generalized time-ordering operator along a closed time contour |^TJ with +(— ) for boson 
(fermion) operators. For the x an d U fields, we have subtracted the classical expectation 
values in order to separate the quantum fluctuations of the fields around the classical field 
configuration. These Green functions generally describe the propagation of an excitation 
in a many-particle system from space-time point x to point y. In the absence of a cou- 
pling between ip, A^ and x, U, and in the zero-density limit, the expectation values Sij 



and can be shown IS] to reduce to the usual quark and gluon Feynman propagators, 
respectively, and A with the scalar Feynman propagator. 

Using the definitions fl27p and implementing the gauge-fixing constraint in (^) for the 
gluon fields, one finds from ( f25| ) the following equations of motion for the Green functions 
(in the limit of zero rest masses), 



ij ■ d x Sij(x,y) = 6ij5 (x,y) + / d x E ik (x, x) S kj (x', y) 

d 2 x D%(x,y) = 5 ab 6\x,y) (V - - |dVll£ a>b ,(x,x') D^ b (x',y) 

d 2 x A(x,y) = -5\x,y) + J d 4 x' E(x, x') A(x', y) 

d 2 x A(x,y) = -8\x,y) + J d A x' E(x, x') A(x' , y) , (28) 

describing the change with respect to x, plus similar equations for the change with y by the 
substitions d x — > —d y , and S(x, x')S(x', y) — > S(x, x')H(x', y), etc.. The explicit expressions 



for the self energies S, II, H, and 3 are rather lengthy and can be found in Ref. [19]. We 
remark that the functions E and II include both the usual quark and gluon self energies, as 
well as the additional interaction of the quanta with the x field. Similarly, the self energy S 
incorporates the effective self interaction of the x field described by the potential ([ll]), plus 
the interaction with the quark and gluon fields in (^). Finally, the self energy of the U 
field plus its coupling to the x field is described by 3. In the equation for the gluon Green 
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function D^ b , we have on the right-hand side the remnant of the gauge constraint, 



£, 



[IV 



d A p e~ ip '( x ~ y } 



(2vr) 4 p 2 + i0+ 



s=l,2 



(29) 



which involves a sum over the two physical (transverse) gluon polarizations (e.g. in Feynman 



gauge e 



A - 5 * and thus £ s 



0). We note that equations (28) are 



of the form of Dyson- Schwinger equations [22], and can be rewritten in symbolic operator 
notations as 



S — Sq + Sq S S , 
A = A + A SA , 



D = D - D UD 
A = A + A HA 



(30) 



where So, Do, Ao, Ao denote the free-field Green functions that satisfy the equations of 
motion in the absence of self and mutual interactions. Fig. 3 illustrates the diagrammatic 
representation of the Green functions S, D, A, A, the self-energies E, II, H, H, and the Dyson- 
Schwinger equations (|3(i|). 

A quantum transport formalism can be derived from the equations ( f28| ) that is very 
suitable for the present purposes |jT9f| . We confine ourselves here to sketching the essential 



steps. One introduces the Wigner transforms W [23] of the Green functions and the self 
energies W = S, D, A, S, II, S: 



W(r,p) 



<?Re ip R W{r, R) 



(31) 



where 

W(r,R) = w(r + ^,r-^\ = W(x,y) , (32) 

with r = {x + y)/2 and R = x — y denoting the center-of-mass and relative coordinates, 
respectively, and R being the canonical conjugate to the momentum p (as before r, p, etc., 
denote four vectors, and a ■ b = a^b^). The equations of motion for the Wigner transforms 
W(r, p) are now obtained |n| under the assumption that the Green functions and self- 
energies W(r, R) can be approximated by a gradient expansion in r up to first order: 



d 

W(r + R,R) ~ W(r,R) + R ■ —W(r,R) . 

or 



(33) 



This assumption implies a restriction to quasi-homogenous or moderately inhomogenous 
systems, such that the Green functions vary only slowly with r. In homogenous systems, 
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such as the vacuum, translation invariance dictates that the dependence on r drops out 
entirely, and the Wigner transforms then coincide with the momentum-space Fourier trans- 
forms of the Green functions and self energies. Although we will in the present paper 
consider the evolution of a fragmenting jet system in vacuum, the subsequent formulation 
is tailored to apply also to more general translation non-invariant situations in moderately 
inhomogenous media. 

The Dyson-Schwinger equations ( f28| ) can now be converted into kinetic equations by 
performing the Wigner transformation for all Green functions and self energies, and 
using (|33|). One arrives at two distinct equations for each of the Wigner transforms S, T> 
and A (A), with rather different physical interpretations: (i) a transport equation, and (ii) 
a constraint equation. The transport equations are 



'■ - •■ 1 — ~(+) i a *•(-) 



p-d r S t ,(r,p) = -Sij-y-dr + i ( 7 .p + E)jrW + - 7 ■ d r ^ 

P-d r v%(r, P ) = ig£<+> 



p ■ d r A(r,p) 



1 



(34) 



whilst the constraint equations are 
1 



P 



S kj (r,p) = <% ( 7 • p + E) + i( 7 - p + E) J* > 



1 



- -X) Sac ~ K,a,c(r,P) 



V 



-d 2 r + E(r,p) 



A(r,p) 



1 



- -H-. 



>ab 



+ 7 Gab (_) 



(35; 



The equations for A are formally identical to those of A. The operator functions Q, 
and 7i (T~i), which include the effects of spatial inhomogenities, are given by (df? = d/dr^ 1 
d£ = d/dp"): 



{-) 



r uu (-) 

y a b 



H~ 



d^E,d r A 



d?Z,d?S 



ab 



•'ij 



{d^,d;s} + -{d^,df l s} 
= ({d»n,d;v} + -{d?n,df l v} 



y ab 



Eqs. (|34j ) and (|35| ) are our general master equations. The physical significance QlSj] of 
the transport equations (|34|) and the constraint equations (|35|) is that the former essentially 
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describe the space-time evolution of the Wigner transforms, whereas the constraint equa- 
tions describe the "orthogonal" evolution in momentum space, and express a normalization 
condition imposed by unitarity and the renormalization group. In order to relate these 
operator equations to physically-relevant (observable) quantities, we define the Wigner op- 
erators F a (r,p) (a = q, g, x, U) in terms of the operators S, V, A, A and the self energies 
S, IT, A, A as follows: 

i5ij(r,p) = 5i d (7 -V + E) (2tt») 5 (f - Y 2 ) F q (r,p) 

W%(r,p) = 6 ab e^(p, S y a *(p,s) (2 m )s{p 2 -Il) F g (r,p) 

*A(r,p) = (2m) 5 (p 2 - a) F x (r,p) 

»A(r,p) = (2tt») <5 {f - 3) Fy(r,p). (37) 

Then, by tracing over color and spin polarizations, and taking the expectation values (or, 
in a medium, the ensemble average) of these Wigner operators, one obtains the scalar 
functions 

F a (r,p) = F a (t, f;p,p 2 = M 2 ) (a = q, g, X , U) (38) 

with 

F q (r,p) = (Tr[S(r,p)}) , M 2 = S 2 (r,p) 

F g (r,p) = (Tr[V(r,p)]) , M 2 = IL{r,p) 

F x {r,p) = (A(r,p) ) , M 2 = E(r,p) 

Fu(r,p) = (Tr[A(r,p)}) , Mjj = S(r,p) . (39) 

The c- number functions F a (r,p) are the quantum- mechanical analogues of the classical 
phase-space distributions that measure the number of particles at time t in a 7-dimensional 
phase-space element d 3 rd p. Due to the effects of the self energies, three-momentum and 
energy are generally independent variables, because the quanta can be off mass shell, i.e., for 
zero rest masses, E 2 = p 2 + M 2 ^ p 2 , where M 2 , eq. (|39|) , represents the off-shellness due 
to the self and mutual interactions of the quanta (M 2 = for on-shell particles). In contrast 
to the classical propagation of on-shell particles, the Wigner functions ( |37| ) incorporate the 
quantum "Zitterbewegung" even in the absence of interactions with other particles. These 
spatial fluctuations arise from the combination p 2 — d 2 /4 acting on the Wigner operators in 
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( p5D , and account for the uncertainty in spatial localization of a quantum particle due to its 
momentum that is determined by the particle's self interaction. The self-consistent solution 



of the transport and constraint equations (34) and (35) therefore corresponds to summing 



over all possible quantum paths r in space-time with fluctuations in energy-momentum p, 
constrained by the uncertainty principle . This simultanous evolution in r and p of the 
Wigner functions F a is illustrated in Fig. 4. 

3.3 Macroscopic quantities 

The functions F a , eq. (^), contain the microscopic information that is required for 
the statistical description of the time evolution of a many-particle system in complete 



phase space [24[. Depending on the physical situation under consideration, one starts 
from specified initial distributions at t = to and follows the time evolution of the phase- 
space densities F a (r,p) according to the master equations (|34]) and (|35|). At any time 
t > to, F a (r,p) reflects the state of the system around r and p. One can then calculate, 
in a Lorentz-invariant manner, directly from the microscopic densities F a , the relevant 
macroscopic quantities that are related to observables. Relativistic transport theory |p5| 
relates physical quantities to phase-space integrals over products of combinations of four- 
momenta or tensors and the particle distributions. Specifically, using (|37|), performing the 
traces over color and spin indices where necessary, and taking the ensemble average, one 
obtains the local space-time-dependent particle currents n a and the corresponding energy- 
momentum tensors Tt u for the different particle species a = q,q,g,X>U, which are given 



by m 



<(r) = J dn a p»F a (p,r) , 2f(r) = J dn a p»p u F a (r,p) , (40) 

where dVL a = r y a dM 2 d 3 p/(l6ir 3 p ), the 7 Q are degeneracy factors for the internal degrees 
of freedom (color, spin, etc.), M a measures the amount by which a particle a is off mass 
shell result of the self energy terms in fj§), @, and p° = E = + v / p 2 + M£. These 
macroscopic quantities can be written in Lorentz-invariant form by introducing for each 
species a the associated matter flow velocity u{J(r), defined as a unit-norm time-like vector 
at each space-time point, (u^u^) a = 1. A natural choice is e.g. = n^/^/n ua n'^. By 
contracting the quantities (|40| ) with the local flow velocities it£J, one can now obtain corre- 
sponding invariant scalars of particle density, pressure, and energy density, for each particle 
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species a individually. For instance, in the absence of viscosity: 



n a (r) 



P a (r) 




-|w) (V- 

T^air) <(r)<(r) 



- u> 




(41) 



Due to the scalar character of these quantities, they provide local, Lorentz-invariant mea- 
sures of the many-particle system. It is left to convenience in which Lorentz frame the 
calculation is performed, but in general it is considerably simplified in the local rest frame 



a given time can then be obtained by integrating over position space. 

4. PARTON-HADRON CONVERSION OF FRAGMENTING JET SYSTEMS 

The preceding kinetic formulation allows us now to apply the conceptual ideas of the 
effective field theory of Sec. 2 to the dynamics of parton-hadron conversion in rather general 
situations. In accord with the formalism of Sec. 2, the parton-hadron transition can be 
visualized as the conversion of high-momentum colored quanta of the fundamental quark 
and gluon fields into color-neutral composite states that correspond to local excitations of 
the condensate fields x an d U embedded in the physical vacuum. 

Ultimately, we would like to address the dynamics of the (non-equilibrium) QCD phase 
transition in finite-temperature systems. Here, however, we will as a first application 
study a much simpler sytem, namely the fragmentation of a qq jet system with its emitted 
bremsstrahlung gluons, and describe the evolution of the system as it converts from the 
parton phase to the hadronic phase (illustrated in Fig. 5). A time-like virtual photon 
in an e + e~ annihilation event with large invariant mass Q ^> Aqcd {^-qcd = 200 — 400 
MeV), corresponding to a very small initial size L <C L c (L c = 0.5 — 1 fm), is assumed to 
produce a qq pair which initiates a cascade of sequential gluon emissions. The early stage 
is characterized by emission of "hot" gluons far off mass shell in the perturbative vacuum. 
Subsequent gluon branchings yield "cooler" gluons with successively smaller virtualities, 
until their mutual separation approaches L f» L x . As is evident from the previous Fig. 2, 
this point characterizes the beginning of the transition, because the partons can now tun- 
nel through the developing potential barrier and form color-singlet composite states, which 



of the matter where 



(1,0). The total number and the free energy of the particles at 
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represent the particle excitations of the scalar long-range x held. These "pre-hadronic" ex- 
citations must then convert into physical hadronic states — either excited gluonic states, or, 
via coupling to the U field, quark-antiquark meson excitations — and subsequently decay 
into low-mass hadrons. 

4.1 General concept 

As stressed in Sec. 2, the phenomenological color-singlet function % represents the effect 
of the long-range order of the non-perturbative vacuum, so that the formation of a color- 
neutral parton cluster - or bubble - can be interpreted as a domain structure immersed 
in the medium of the non-perturbative vacuum. In quantum field theory such stable field 
configurations arise as classical soliton solutions of the equations of motion || 26 ] . On the 



other hand, it is well known that QCD exhibits the so-called 'preconfinement' property p^ ] 
already on the perturbative level, which is the tendency of the gluons and quarks produced 
in parton cascades to arrange themselves in color-singlet clusters with limited extension in 
phase space. It is therefore natural to suppose that these clusters, or bubbles, are the basic 
'pre-hadronic' units out of which hadrons arise non-perturbatively. 

Thus, the kinetic evolution of the system develops in three stages: parton multiplication, 
parton-cluster conversion, and cluster decay into hadrons. It is clear that in this approach 
the conversion process is a local, microscopic mechansim, that proceeds earlier or later 
at different points in space, depending on the local density of partons and their nearest- 
neighbour separation L, as defined by eq. (|3|). Thus, in order to trace the full dynamics, it 
is necessary to follow the evolution of the particle distributions in real time using the kinetic 
framework of Sec. 3. In accord with the above picture, we will now proceed in several steps, 
starting from the master equations (Q) and p5|): (i) employ a separation of a coherent 
(mean) field part and a contribution from quantum excitations for the composite fields x 
and U, (ii) fix a specific, ghost-free gauge for the gluon fields that is most convenient for our 
purposes, (iii) treat the evolution of the high-momentum quarks and gluons perturbatively 
in the presence of the coherent field x- 

(i) According to our interpretation of oscillations about the minimum of the potential at 
Xo, Uq as physical excitations of the coherent fields % an d U, we separate in a standard way 
[p8| the classical field configuration at the minimum of V x=Xo in Fig. 2, from the quantum 
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fluctuations around this minimum. We represent x = X + X an d U = U + U , where 
X, U are c-number functions (the mean field parts), and x, U denote quantum operators 
(describing the excitations). The physics behind this separation is that the coupling of x 
to ip, tp, A^, as well as to U, will make the composite fields x anci U dynamical variables, 
so that the fluctuations around the mean fields x an d U wm propagate and form collective 
excitations. Therefore the system is characterized (aside from the elementary fields ip, 
tp, A^) by the mean fields, as well as by the collective excitations with their own energy 
spectrum and distribution. With this prescription we can treat the local interaction of the 
partonic fluctuations with the coherent field analogously to the familiar problem of quantum 
fields (ip,^,A^) interacting with a classical "external" field (%), which converts the partons 
to color-singlet clusters or bubbles corresponding to excitations in the coherent field (x). 
Specifically, in our approach the bubbles represent non-topological soliton configurations 
which are stable, classical solutions of the equations of motions, as have been studied for 
instance by Friedberg and Lee and Coleman [26]. We do not include here the possible 



additional interactions between partons and bubbles, or among bubbles themselves. 



(ii) It is convenient to work in a physical (axial) gauge p9l 3C] for the gluon fields, generi- 
cally given by choosing the gauge function (, a {A) in (||) as 

UA) = - d ^ n ■ A *) d \ n ■ A a) > ( 42 ) 

where a is the gauge parameter, and is a constant vector with n 2 ^ 0. In particular, 
we will set a = 1 which is known as the planar gauge. In contrast to covariant gauges 
where £ a (A) = —l/(2a)(d ■ A a ) 2 , the class of gauges ( fl2"| ) is well known to have a number 
of advantages. It is ghost-free, i.e. the ghost field contribution in ( |l9| ) decouples and drops 
out. Also, the so-called Gribov ambiguity is not present in this gauge. Another feature of 
( f42| ) is that the gluon propagator involves only the two physical transverse polarizations, 
so that the equations (|34|) and (|35|) simplify considerably [19]. Furthermore, it allows for 



a rigorous resummation of the perturbative series at high energies in terms of the leading 
logarithmic contributions and consequently leads to a simple probabilistic description of the 
perturbative parton evolution within the (Modified) Leading Log approximation (MLLA) 



[]31|, |3| in QCD. 

(iii) We will evaluate iteratively the 2-point Green functions S and T> of quarks and gluons, 

22 



respectively, in the one-loop approximation in the framework of "jet calculus" [33], using 
the MLLA techniques of coherent parton evolution including soft-gluon interference [31]. 
The associated quark and gluon self energies S and II include both the one-loop quark- 
gluon self interaction through real and virtual emission and absorption, and the effective 
interaction with the confining background field \ B- Similarly, the self-energy H of the 
X field embodies the self interaction and the coupling to the U field, as contained in the 
effective potential (f[l]), as well as contributions from quark and gluon recombination to \ 
excitations. Correspondingly, the function H of the U field incorporates its self interaction 
and the interaction with the x field. 

4.2 The kinetic equations for real-time evolution in phase space 

As a consequence of the prescriptions (i)-(iii), and of exploiting in the present e + e~ — > 
hadrons case the special property of translation invariance of the parton evolution in the 
perturbative vacuum, one finds after a lengthy calculation that the transport equa- 
tions (|34|) and the constraint equations (|35|) can be combined in a single set of coupled 



integro-differential equations for the phase-space densities F a (r,p) defined by (37 ) and 
(|38|). Introducing the usual light-cone variables 



p V = (j,+ :P rf ± ) 5 P ± = P0 ± Pz , p ± = (p x ,Py) (43) 

and 

x = ^- , p± = \JpI+pI , P 2 = Vfjf ■ (44) 

where Q is the hard scale of the initial qq pair created by the photon, and r = = (t,r), 
we write 

F a = F a {r,p) = F a (t,f;x,pl,p 2 ) . (45) 

The kinetic equations which one obtains from the transport equations ([P|) by implementing 
the constraints ( |35| ) can now be summarized compactly as follows (see Fig. 6): 

tF q = + A q q 9 F q + AfF g - B%F q F g - B% F q F q (46) 

tF g = +AfF g - AfF g + £ A^F q+ - q - B™ F g F g - £ Bf q F g F q+q (47) 

/ / 
KF X = +C^F g F g + Y^Cf q F g F q+q + Y,C^F q F q -t) u x F x -E h x F x (48) 

/ / 

tF v = +D' U X F X - E>%F V , (49) 
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where we have abbreviated F q+q = F q + F q , which in the present case is equal to 2F q , 
and the / means summing over quark flavors / = u,d,s, . . .. The left-hand sides of these 
equations describe the propagation of the particles in the presence of the mean field, whereas 
the terms on the right-hand sides represent the effects of particle creation, annihilation, and 
recombination. The operator fC on the left-hand sides is given by 



ICF a 



p^d? + (M a d?M a )dP 



F a , (50) 



with the first term describing the free propagation and the second term reflecting the effect 
of the mean field (Fig. 6a). The functions M a are the mean-field parts of the self energies 
M a defined in (|39|). On the right-hand side the quantities I^i'^ 2 'Z (where I = A,B,..) 
represent integral operators that incorporate the effects of the self energies in terms of 
the relevant amplitudes for the various processes a±,a2, ■■ — ► &i,&2> ••> and that act on the 
phase-space densities to their right (Fig. 6b). These coupled equations reflect a probabilistic 
interpretation of the evolution in terms of successive branching and recombination processes, 
in which the changes of the particle distributions on the left-hand sides are governed by the 
balance of gain (+) and loss (— ) terms on the right-hand sides. The different terms on the 
right-hand sides of eqs. (^)-(]49|), the contributions to the gain and loss of particles, fall 
into three categories: 

(a) Parton multiplication through emission processes q — > qg, g — > gg and g — ► qq; 

(b) Parton cluster formation through recombinations qq — ► XX> Q9 ~> QXi 99 ~^ XX] 

(c) Hadronic cluster decay either through direct conversion of the formed scalar \ excita- 
tions into hadrons h through x — ► h, or via coupling to the pseudoscalar states \ ~~ > U, and 
the subsequent decay into hadrons, U — > h. 

In the following subsections we explain these contributions in detail. 
4.3 Parton multiplication 

The integral operators A in the quark and gluon evolution equations (|46|) and 



represent the changes of the parton distributions in phase space due to the perturbative 
cascade evolution. Explicitly, 



AfF q = X x I dz 



Jo dZ \ F ^{^'^ zp2 ^ zp2 ^) ~ F g( r '^ x 'P±'P 2 ) 
AfFg = X x y F g (r]^,zp\,zp*\ -y g ^ qq -(z,e) a g {z,p 2 ) 



lq^qg{z,e) CL q (z,p 2 ) 
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A qq F 



dz 



F g r; -,zp ± ,zp' 



- F g (r;x,p ± ,p ) 



X x n f (p ) F g (r;x,p ± ,p ) / dzj g -, q g(z,e) a g (z,p ) 

•A) 



A q q F q+q = X X J Q ~ F q+q [r\ ~, Zp ± , Zp ) J q ^ gq (z , e) (l q (z , p ) . (51) 

Here A x = A x (x(r)) = 1 - (x/xo) 4 + O[(x/xo) 6 ], and the function a(z,p 2 ) is given by 

a g , s (^,p 2 ) := i- T g , g (p 2 ) a, ( (1 - z)/) , (52) 

with a "life-time" factor T q ^ g {p 2 ) that expresses the probability for a parton of virtuality p 2 
to decay (branch) within a time interval i in the laboratory frame, 

t 



TUP" 



1 



cxp 



T qAP 



where r{p 2 ) oc E/p 2 (explicit expressions can be found in Ref. |3 
the one-loop QCD coupling 
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a s (k 



33-2n f (k 2 )) ln[(k 2 + k 2 )L 2 ] 



(53) 

Furthermore, a s is 
(54) 



and nf{k 2 ) is the effective number of quark flavors at k 2 



n f (k 2 



/ v 



Am 2 



(55) 



In (|54|) we have assumed the correspondence L c ~ A.qq D to the intrinsic perturbative QCD 
scale, and k$ is a parameter that prevents a divergence when k 2 — ► £,T 2 , and defines a 
maximum value a s (0). We will determine ko in Sec. 5 from the total parton multiplicity. 
The functions 7(2, e) are are analogous to the standard branching kernels in the MLLA |29[ |. 
Note that the 4-gluon vertex vertex does not contribute in the MLLA in the gauge 
because it is kinematically suppressed. As a consequence, the effect of the couplings kl(x) 
and n(x), eqs. ([l(]) and (17), on the parton evolution reduces to 2-parton recombinations 
into color-singlet clusters — the terms proportional to B, C which will be given below. 

The branching kernels ^f a ^ c {z) are the familiar energy distributions for the branching 
a — > be with z = Xb/x a and 1 — z = x c /x a the energy fractions of of daughter partons: 



-y q ^ qg (z,e) 



C F 



1 + z 2 
I- z + e 
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1g-+ggM = 2C A (z + ^—-^ + z{l - , 

\1 — z + e z + e 

7 9 ^,e) = i (z 2 + (l-z) 2 ) , (56) 

where Cp = (N% — 1)/(2N C ) = 4/3, Ca = N c = 3. In the denominator of 7 g ^ g9 and 7 ff ^ g9 , 
there appears the function 

'22 2 

4(p • n) 2 p2 v ) 

where p (p') is the momentum of the mother (daughter) parton and p± the relative trans- 
verse momentum of the daughter partons with respect to the mother. It arises here as a 
consequence of the constraint equations (|35]) which determine spatial uncertainty associ- 
ated with the off-shellness of the partons. It effectively cuts off small- angle gluon emission 
by modifying the free gluon propagator oc Zg 1 to the form (z g + e) _1 (where z g = z or 
z g = 1 — z) when p±/p z = 0(1), that is, in branching processes with large space-time un- 
certainty. This ensures that the two daughter partons can be resolved as individual quanta 
only if they are separated sufficiently by Ar± oc 1 /p± in position space, in accord with the 
uncertainty principle. Note that e can be neglected in the terms oc (z g + e) _1 in (|56| ) for 
energetic gluon emission (z g — * 1), but is essential in the soft regime (z g — > 0). The effect of 



e has been shown [29, |35|] to result in a natural regularization of the infra-red-divergent be- 
haviour of the branching kernels (^) , due to destructive gluon interference which becomes 
complete in the limit z g — > 0. 

4.4 Parton cluster (bubble) formation 

The operators B, C in eqs. fl46|)-(|4"8|) represent the changes of the phase-space densities 
due to recombinations of two partons at r and r' to color-neutral clusters, or bubbles 
that arise as non-trivial structures in the vacuum because of the confinement mechanism. 
Their formation is determined, in analogy to the finite-temperature QCD phase transition 



[14|, by the probability for tunnelling through the potential barrier of V between x = an d 
X = Xc in Fig- 2, which separates perturbative and non-perturbative vacua. The associated 
rate of bubble formation around L = L c is generically given by an exponential probability 
distribution (36|, [37]], 

tt(L) = ML) ( 1 - exp[-AF L]) , (58) 
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where ttq = const. ln(l — u)9(l — 2u) + 9(2u — 1), with u = AF L, modifies the small-L 
behaviour for which the exponential form is not appropriate. Here AF is the change in the 
free energy of the system that is associated with the conversion from partons to clusters. 
In our case (for baryon-free matter in general), 

AF = E vol + E surf = — R 3 (L) AP(L) + 4vr R 2 {L) a(L) , (59) 

where R(L) is the radius of the bubble depending on the parton separation. The first term 
is the volume energy determined by the difference of pressure in the perturbative and the 
non-perturbative vacuum, 

AP(L) =P qg (L) -P X (L). (60) 
The second term in (|59|) is the surface energy with the surface tension estimated to be 



a(L) = r d X y/2V(L) « 2 d X ^2V(L) (61) 

"0 J Xmax 

where Xc corresponds to the local minimum of V at L c and Xmax is the point of the local 
maximum of V that separates unconfined and confined domains, as defined in Fig. 2. 

A parton bubble is stable if dAF/dR = 0, which leads to the condition for the stationary 
bubble radius, 

(62) 



R c - R(L C ) = |^ 



L=L C 



When inserted in (|59D, this gives for (pq) 

47T 



7r(L) = 7r (L) 1 - exp 



Rla c L 



3 -^c 



(63) 



with a c = a(L c ). In accord with our definition (^), we interpret the space-time scale L 
as the characteristic inter-parton separation, that is, we define it in terms of the distance 
measure Ajj between two partons, labeled with indices a and 6, 



A <*6 = V r ab r ab,^ , r ab = r a - r b , (64) 

and identify L with the the minimum distance L ab for a certain parton a to its next 
neighbour b: 

L(r) = L ab = min 6 (A al , . . . , A a6 , . . . , A an ) . (65) 

Other measures are possible, e.g. A ab = (l/x^ b + + + 1/(^) 2 ) _1 ^ 2 j but we find 
that the particular choice of A ab is not crucial as long it provides a reasonable distance 
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measure. We choose (|64|), because it has the advantage that it is manifestly Lorentz- 
invariant and has a simple interpretation as the two partons' spatial separation in their 
center-of-mass frame. 

We assume that the dominant contribution to bubble formation arises from 2-parton 
fusion and ignore recombinations of 3 or more partons. This is reasonable, unless the local 
parton density becomes so large that also the latter processes have a significant proba- 
bility to occur. We take the same probability distribution tt(L) for the various types of 
configurations, since it depends only on the 'color- and flavor-blind' variable L, i.e. we set 

n(r,r') = vr(L) = ir gg ^ xx (L) = vr 9 ^ xx (L) = ir gq -> xq (L) , (66) 

where L is given in terms of r and r' by (|65|). The various B terms in (^) and (^7|) can 
then be expressed generically as: 

B%F a F b = F a (r;x,pl,p 2 ) J dVvr(r,r') J d*p' F b (r';x',p'?,p' 2 ) (67) 

where a,b = q,q, g and c = x, g, and d 4 p' = dp 2 dp 2 dx' /x' . Similarly the C terms in fl48[) 
are given by: 

C%F a F h = Jd 3 r' Jd 3 r"Tr(r>,r")5 3 (r-^ 



4.5 Hadronic cluster decay 

The ensemble of clusters determined by the coupled equations (|46|)-(|4S|) yields a conti- 
nous mass spectrum of color-singlet excitations with different flavor contents corresponding 
to the types of recombined partons. These states must then decay into physical hadronic 
states with a discrete mass spectrum. The invariant mass distribution of the formed clus- 
ters may be interpreted as a 'smeared out' version of the spectrum of primordial resonances 



formed in the early stages of the confinement mechanism |31]. It therefore seems reasonable 
to treat the fragmentation of these central clusters as a kind of averaged resonance decay 
which, as implied by our locality assumption, must be determined entirely by their invariant 
masses, flavors and total angular momenta. Each cluster in the resonance spectrum may 
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either represent a single hadron resonance that converts directly into a physical hadron 
with a definite mass, or else fragments through a two-body decay into a pair of final-state 
hadrons. From the particle spectra obtained in e + e~ -annihilation experiments it appears 
that quasi-two-body final states are universally dominant, so that the latter possibility 
seems favored if kinematically allowed. 



We adopt the cluster fragmentation scheme of Refs. [38, [39[, however with some mod- 
ification concerning heavy clusters. We assume that each cluster C = x>U can decay by 
either one of the following mechanisms: 

(i) If a cluster C is too light to decay into two hadrons, it is taken to represent the 
lightest single hadron (meson) h, corresponding to its partonic constituents, C — > h, with 
its mass shifted to the appropriate value by adjusting its energy through exchange with a 
neighbouring cluster. 

(ii) If, however, a cluster is massive enough to decay, it decays isotropically in its rest frame 
into a pair of hadrons (mesons or baryons), C — ► h± + h,2 according to the decay probability 
specified below. 

Occasionally it occurs that a cluster comes out very heavy, in which case isotropic 2-body 
decay is not a reasonable mechanism any more. In this case we impose the constraint that, 
if a cluster is heavier than a critical threshold M cr a = 4 GeV, then it is rejected and the 
two recombining partons of that potential cluster propagate on as individual quanta, and 
continue to participate in the parton cascade process, either until they have decreased their 
virtuality sufficiently, or until they recombine with a lower-mass partner. 

To implement this scheme, we observe from eqs. ( |48D and fl49| ) that the cluster-hadron 
transformation can proceed through the scalar channel x ~ * h, or via the pseudoscalar 
channel x ~> U ~* ^i^2> depending on the corresponding density of states with masses 
below the decaying cluster. We assume a Hagedorn [ 40 1 density of states 



p h (m) = cm a exp ( _ ^rj > ( 69 ) 

where c, a are constants and To is the Hagedorn temperature with the typical values c = 
8m^, a = 3 and To = m^. The decay probability of a cluster with mass mc = \p$? to 
decay into a hadron state of mass m' h = \fp 2 is then given by 

1 o r m c=Vi^ 



^C->h(p,p') = jVc( p 2) T c(P 2 ) I . r^dmphim). (70) 
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Here 

N C (P 2 ) = I ' dm p h (m) , (71) 

J rn.Tr 

and in analogy to (|53"D, Tc(p 2 ) is a "life-time" factor giving the probability that a cluster 
of mass tUq = p 2 decays within a time interval t in the laboratory frame, 

T c {p 2 ) = 1 - exp [ , (72) 

where in this case we simply take Tc(p 2 ) = Eq/p 2 = j/mc from the uncertainty principle. 
In order to find the value for the decay probability ( |70|) for a given cluster of mass mc, 
we sum over the possible decays for this cluster according to the particle data tables. The 
probability for a specific 2-body decay mode is taken to be a product of a flavor, a spin and 
a kinematical factor 1381, 



rc-4/n/i 2 (mc;mi,m 2 ) := P m (m c , mi + m 2 ) P s (ji,j2) Pk(m c , m 1 ,m 2 ) , (73) 
where ji 5 2 (^1,2) are the angular momenta (masses) of the two hadrons hi 2- The factor 



is the two-body phase-space suppression function for the decay. The spin factor 

Ps{juh) = (2ji + l) (2j 2 + l) (75) 

takes into account the spin degeneracy with the allowed spins ji and j'2 of the two hadrons. 
The kinematic factor 

^/X(m 2 c ,m 2 ,m 2 2 ) 
P k (mc,mi,m 2 ) = 5 (76) 

is the two-body phase-space factor, where A(a, 6, c) = a 2 + 6 2 + c 2 — 2 (a& + ac + 6c). 

Thus, with the decay probability Tc~>h of ( |70| ) evaluated in this fashion, the terms 
involving the D and i£ operators in the kinetic equations (Hq) and (H9) can be expressed as 



D V X F X = F x (r-x,p\,p 2 ) J dp' 2 T x ^u(p 2 ,p 2 ) 

E\F X = F x (r;x,p 2 ± ,p 2 ) J dp 2 T x _^ h (p 2 ,p 2 ) 

D' X F X = J dp^dp'^F^x^'^p'^T^uip' 2 ^ 2 ) 

E'uFu = Fu(r;x,pl,p 2 ) J dp 2 Tu^ h {p 2 ,p 2 ) (77) 
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4.6 Method of solution by Monte Carlo simulation 



We can now solve the set of evolution equations fl46|)-([49D by means of a real-time 
simulation in full phase space using the computational methods of Ref. [38]. One starts 



from an initial phase-space density of partons, which in the case of a jet-initiating qq pair 
with invariant mass Q is 

F q+q (t = 0,r;x,pl,p 2 ) = 5 3 (r) 5{f - Q 2 ) 6(x - 1) 8(p\) 5(p 2 - Q 2 ) , (78) 

where we choose the qq center-of-mass frame as our reference frame, and we sum over all 
quark flavors / weighted with a factor wj = e 2 /rif(Q 2 )^Jl — 4m 2 /Q 2 6 (q 2 — Am 2 ^j, that 
accounts for the electromagnetic charge and mass threshold of the initial qq pair produced 
by the photon (or Z°). 

The parton shower development is then followed in a cascade simulation (for details see 
|m| , |38|| ) in the center-of-mass frame: The system of particles is evolved in discrete time 
steps, here taken as At = 0.01 fm, in coarse-grained 7-dimensional phase-space with cells 
Ail = A 3 rA 3 pAM 2 . The partons propagate along classical trajectories until they interact, 
i.e., decay (branching) or recombine (cluster formation). Similarly, the formed clusters 
travel along straight lines until they decay into hadrons. The corresponding probabilities 
and time scales of interactions are sampled stochastically from the relevant probability 
distributions according to Sees. 4.3-4.5. At any time t > we can extract the phase-space 
densities (j38[), F a (t, r,p, p 2 ) = dN a (t)/d 3 rd 4 p of the particle species a = q,g,x,U, and 



with these phase-space profiles we can then calculate, using the formulae ( |40| ) and (41), 
the associated local pressure P(r), particle density n(r) and energy density e(r) for each 
species individually, these being the quantities that characterize the macroscopic state of 
the system at t and f. 

With this concept, we can trace the space-time evolution of the parton-hadron conver- 
sion process self-consistently: at each time step, any "hot" off-shell parton is allowed to 
decay into "cooler" partons, with a probability determined by its virtuality and life time. 
Also in each step, every parton and its nearest spatial neighbor are considered as defining 
a fictious space-time bubble with invariant radius L, as defined by (|65|) . By comparing the 
local pressure of partons P qg (t,f,L) with the pressure P x (t,r,L) that such a pre-hadronic 
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bubble would create instead, we obtain the associated value of the conversion probability 
7r(L) which determines the cluster (bubble) formation rate explained in Sec. 4.5. If the 
partons do convert into a cluster, they disappear from a phase-space cell, and instead the 
composite cluster appears at the same space-time point, from which it propagates on. Oth- 
erwise the partons continue in their shower development. The final decay of each formed 
cluster into hadrons is simulated analogously, except that it does not require the comparison 
of pressures, but is determined by kinematics and the available phase space. This cascade 
evolution is followed until all partons have converted, and all clusters have decayed into 
final-state hadrons. 

As an illustrative example, we show in Fig. 7 the time evolution of the particle density 
profiles of partons (Fig. 7a) and clusters (Fig. 7b) for a jet system with invariant mass 
Q = 100 GeV. It is evident how the system evolves in position space with respect to the 
center-of-mass of the two inital partons as a polar wave front (the pictures are symmetric in 
rj_ = ^Jr% + Ty), with the partons gradually converting to clusters. It is interesting that this 
local excitation of the vacuum due to the injection of the highly virtual dijet system, and 
the subsequent evolution, resemble very much the situation of a stone plunged into water, 
with a well-shaped "shock front" expanding isotropically in the center-of-mass frame. 

5. PHENOMENOLOGY 

In this Section we study the observable implications of our approach to parton-hadron 
conversion, and investigate its consistency with standard particle physics phenomenology. 

5.1 Determination of the potential V(x, U) 

We first need to specify the parameters of our approach. Recall that this phenomeno- 
logical input is contained in the effective long-range potential V(x, U), eq. (M), which 
combines with the dynamical contribution 5V(L,x) to the L-dependent potential V(L), eq. 
(pl|). As L varies, V changes its shape, which affects the dynamical evolution of the system, 
and the latter in turn determines the further change of V. Hence, the details of the dy- 
namics are governed by the choice of parameters in V and thus V. The crucial parameters 
are the bag constant B which defines the vacuum pressure V(0) in the short-distance limit 
L — > 0, and XOi the value of the condensate of x i n the long-distance regime. As indicated 
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in Fig. 2, as L increases, the changing form of V(L) is characterized by three distinct length 
scales: L = L x , the point when partons begin to convert, L = L c , when the pressures of 
partons and pre-hadronic clusters equal each other, and L = Lq, when the transition is 
completed. 

We will fix B and xo> which have well-defined physical interpretations, and then de- 
termine L x , L c and Lq. Although the values of B and xo are n °t precisely known, there 
is agreement of various phenomenological determinations about their ranges: one expects 
[14) = (150 - 250) MeV and X o = (50 - 200) MeV. In the following we adopt two 



representative parameter combinations: 

B 1/4 = 230 MeV X o = 200 MeV 

B 1/4 = 180 MeV xo = 100 MeV . (79) 

Then, with the potential V specified, we can determine the values of L x , L c , Lq from the 
Monte Carlo simulation of the evolution of the system as the scale L changes due to the 
particles' diffusion in phase space. The most interesting quantity is L c , the point which is 
characterized by the equality of partonic and hadronic pressures. As explained in Sec. 4.6, 
we can compute the corresponding pressures P qg and P x from the phase-space densities of 
partons and clusters, respectively. In analogy to Ref. [14], we represent (on dimensional 
grounds) 

P qg (r,L) = a qg (r,L)L-* - B 

P x (r,L) = a x (r,L)L- 4 - V(L) , (80) 

and, because V(x,L)\l=l c = V(Xc,L c ) = (c.f. Fig. 2), we have 

~a qg (r,L c ) - a x {r,L c y li ' 



B 



(81) 



The dimensionless functions a qg and a x are obtained from the numerical simulation and 
are shown in Fig. 8 as a function of time for the above two choices of B and xo m the cases 
of qq and gg jet evolution with Q = 10 GeV and Q = 100 GeV. Plotted are the kinetic 



pressures P(t, L) := a(t, L)L~ A (where t ~ \J L 2 + z 2 ) along the "shock front" of the jet 
profile which is seen in the previous Fig. 7. From Figs. 8a and 8b one observes that (i) 
the pressure evolution obviously depends on the type of the two jet-initiating partons: it 
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decreases more slowly for qq pairs of different flavors than for a gg pair, because gluons have 
a larger emission rate and therefore the two leading gluons evaporate their initial energy 
faster; (ii) the crossover point between the pressures P qg and P x is rather insensitive to the 
choice of L c ; (iii) the crossover is shifted away from t = with increased jet energy Q; (iv) 
at L c the partonic pressure P qg still exceeds P x , i.e. a qg (L c ) > a x (L c ), consistent with (|8l"|). 
From this analysis, we find using the determining condition (|Sl"|), 



0.6 fm forB 1 / 4 = 230 MeV 
L c = { . (82) 

0.8 fm for B 1 / 4 = 180 MeV 



Directly associated with the scale L c is the parton-cluster conversion probability 
which is determined by the width of the potential wall between the two phases. It enters 
the kinetic equations via ( |67| ) and determines locally the time scale of the parton-to-cluster 
transition by the magnitude of the surface tension a c as given by (|6l|). We find 

1/3 J 40 MeV for L c = 0.6 fm {B 1 ' 4 = 230 MeV) 
<* J = \ ... > (83) 

! 48 MeV for L c = 0.8 fm {B 1 / 4 = 180 MeV) 

which by virtue of (63) fixes the cluster-formation probability tt(L). It is noteworthy that 
the above small values of the surface tension a c are in agreement with lattice QCD simu- 



lations [41]. and correspond to a weakly first-order transition at finite temperature, which 
is consistent with astrophysical constraints [42] on inhomogenities. This finding implies a 
rather rapid conversion of partons into color-singlet clusters (pre-hadrons), as is also evident 
from Fig. 8. This means that parton-hadron conversion is not dependent on the details 
of the interpolation functions hl(x) an d fJ>i(x)i as already advertised in Sec. 2.3, and has 
interesting consequences for the cluster size and mass distribution, as we will discuss below. 

The value of L x , below which size only the perturbative vacuum of the pure parton 
phase can exist, is given by the point of inflection of the effective potential V, when the 
local minimum at (x) ^ ceases to exist (c.f. Fig. 2). It turns out to be rather close to L c , 

I 0.4 - 0.5 fm for L c = 0.6 fm (B 1 / 4 = 230 MeV) 

L x ~ \ i ( 84 ) 

I 0.6 - 0.7 fm for L c = 0.8 fm (B 1 / 4 = 180 MeV) 

which is a consequence of the small values ( |83| ) for the surface tension a c , and indicates that 
the transition occurs very abruptly. Finally, we find that the scale Lq, when the parton- 
cluster conversion is complete and no partons are left over, depends not only on L c but also 
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on the initial jet energy Q. It gives an estimate for the time scale tq oc Lq of the global 
conversion process, and comes out rather long, namely for L c = 0.6 fm we get To ~ 8.5 
(21) fm for Q = 10 (100) GeV, whilst for L c = 0.8 fm the time scale is t « 10 (26) fm 
for Q = 10 (100) GeV. 

Immediate consequences of the values for i? and xo m ( f^9[ ) are the "critical temperature" 
for the phase transition in finite-temperature QCD, 



9B\ 1 / 4 I 160 MeV for L c = 0.6 fm (B 1 / 4 = 230 MeV) 
47f2 / I, 125 MeV for L c = 0.8 /m {B l > 4 = 180 MeU) 

the characteristic mass scale of the lightest scalar glueball, given by [0] 



(85) 



ld*V( X ,Q) 




dx 2 


X=X0 



i ^/B I 1.05 GeV for L c = 0.6 fm 
I 1.30 GeV forL c = 0.8 fm 



and, the estimate for the value of the gluon condensate (23) 



32 1.25 GeV fm" 3 for L c = 0.6 fm 

Go = 4 B = ■ ( 8? ) 

y I 0.50 GeV fm~ 3 for L c = 0.8 fm 

The parameter values obtained above are summarized in Table 1. Both choices (|7^) of B 1 ^ 
and xo give reasonable results that are in the range of commonly-accepted phenomenology. 



5.2 Cluster distributions and hadron spectra 

Using the parametrizations of Table 1, we have investigated more quantitatively a num- 
ber of typical features of the jet evolution, which we discuss now. 

In Fig. 9 we show the total transverse momentum generated during the time evolution 
of the system in the center-of-mass of the initial jet pair: 

p W(t) = Jd 3 rJ dxd P 2 dp 2 ±P± F a (t,f;x,p 2 ± ,p 2 ) , (88) 

where a labels 'partons' or 'clusters', and p± = ^Jp 2 + p 2 . As before, we compare the cases 
Q = 10 (100) GeV and L c = 0.6 (0.8) fm. At t = we start with p^\o) = 0, because the 
two initial partons recede back-to-back along the z axis. Then, with progressing time, the 
jet evolution can roughly be divided into four stages: (i) a very short hard stage ( <■ 0.02 
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fm) , characterized by an explosive production of partons and consequently transverse mo- 
mentum; (ii) a longer shower stage (~ 0.02 — 0.3 fm) that essentially just causes diffusion 
in phase space with little additional entropy and transverse momentum production; (hi) a 
conversion stage (~ 0.3 — 5 fm) which sets in when the partons start locally to form clus- 
ters; (iv) a hadronization stage ( > 5 fm) when the clusters start fragmenting into physical 
hadron states via cluster decay and resonances. 

It is interesting to inspect the distribution of the cluster sizes and the invariant-mass 
spectrum of clusters, since these measures are essentially the only microscopic informa- 
tion (aside from the momentum spectrum) which is carried over from the partonic to the 
hadronic phase, and therefore determines directly the final-state hadron distributions. In 
Figs. 10 and 11 we show these distributions, for the two values of L c and the two jet 
energies Q considered before. 

From Fig. 10a it is obvious that the typical cluster radius is strongly peaked at a value 
slightly above L c , with a very small width of about 10 _1 fm. Specifically, the average 
cluster size is {R c i} = 0.62 (0.83) fm for L c = 0.6 (0.8) fm. The narrow width of the 
cluster size distribution is a consequence of the smallness of the surface tension a c (|83|), 
which implies a very small surface energy to be overcome, and results in a jump in the 
tunnelling probability (^) around L = L c . One also sees that there are a few clusters with 
radius smaller than L c , a feature which arises from the fact that cluster formation sets in 
already at L x (c.f. Table 1), when there is a non-zero transition probability for partons 
to yield clusters by tunnelling through the potential barrier that begins to develop when 
L > L x (c.f. Fig. 2). A striking feature of the cluster mass spectrum shown in Fig. 10b 
is that it is insensitive to the choice of L c . The spectrum is characterized by a strongly- 
damped exponential form at low mass (M c i <, 3 GeV), with a high- mass tail that extends 
substantially beyond 5 GeV. For both choices of L c , the value for the average cluster mass 
is (M c i) ~ 1.2 GeV, which is in qualitative agreement with the characteristic mass scale 
m x obtained in (jjffi ) as an estimate. 

Figs. 11a and lib display the corresponding size and mass distribution for the case when 
we impose the constraint mentioned in Sec. 4.3, that two partons cannot form a cluster if 
their invariant mass is above a maximum cluster mass M cr n, even if their mutual separation 
increases beyond L c , where we chose M cr n = 4 GeV. In this case the two partons propagate 



36 



on as individual quanta until they have either radiated off further energy, or recombine with 
a less-energetic partner. Consequently, the cluster size distribution is shifted to larger radii, 
not however very significantly, because only a small fraction (typically less than 5 %) of the 
clusters are affected. At the same time, the mass distribution loses its high-mass tail and 
falls off considerably above 3 GeV, and the average cluster mass comes down substantially 
to (M d ) ~ 0.8 GeV. 

The most remarkable result of Figs. 10 and 11, however, is that the shapes of both the 
cluster size and cluster mass distribution are essentially independent of the jet energy, as 
well as of the initial 2-jet configuration, and therefore appear to be universal. 

The decay of the spectrum of formed clusters into hadrons, simulated according to 
the Sec. 4.5, then yields the average particle multiplicities of final-state hadrons. It is 
interesting to look at the relation between parton and hadron multiplicities. The feature 
evident in Figs 10b and lib, namely that the mass spectrum of color-singlet clusters is 
independent of the total jet energy, is in agreement with analytical predictions |3l}| . This 
implies that parton and hadron multiplicities should be proportional to each other at high 
energies, which is known as local parton-hadron duality |32[| . Fig. 12 display our results 
for the total gluon and quark multiplcities {n qg ) = (n g ) + J2f( n q + n q)-> as wen as the 
ratios of charged hadrons to partons (n c h) / {n qg ) , and of clusters to partons, (n c /) / (n qg ) , as 
a function of jet energy Q. In Fig. 12a the calculated rise of (n qg ) is shown for the two 
choices of L c . The smaller value of L c gives a larger multiplicity because we identified L^ 1 
with the scale Aqcd in ol s , eq. (|54"|). From Fig. 12b one reads off, however, that the average 
number of clusters per parton is independent of L c and of energy Q. As a consequence, the 
number of charged hadrons per parton is larger for larger L c , because at any fixed Q and 
due to 4-momentum conservation, fewer partons yield more massive clusters which in turn 
decay in a larger number of low-mass hadrons. Most important, one sees that the ratio 
i n ch} I ( n gg) is onr y weakly energy-dependent for Q > 30 GeV, in accord with the hypothesis 
of local parton-hadron duality: it rises over this range by less than 10% and appears to 
saturate asymptotically, approaching a constant of ~ 1.6 (1.7) for L c = 0.6 (0.8) fm. 

The resulting average multiplicity of charged hadrons as it rises with Q is shown in 



Fig. 13a together with experimental data [45, 46]. In order to obtain the correct overall 



normalization, we fitted the the infrared regulator ko entering a s , eq. (54), to give the 
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measured total charged multiplicity at Q = 91 GeV for L c = 0.8 fm. The required value 
is ko = 0.5 GeV which implies a s (0) ~ 0.8. With this adjustment we then obtain in Fig. 
13b the momentum spectra of charged hadrons with respect to the variable ln(l/x), where 
x = 2E/Q is the particle energy normalized to the total energy Q. The spectra clearly 
exhibit the well known 'hump-back plateau' |3^]. The good agreement of the simulation 
with experimental spectra is another indicator of the aforementioned local parton-hadron 
duality. The result is not surprising, since the simulation incorporates the coherent parton 
shower evolution ^] based on the soft-gluon interference properties of the MLLA, which 
cause this hump-back plateau with its depletion at small x. Because in our approach 
cluster formation and subsequent cluster decay involve no momentum dependence, but are 
solely described by the space-time separation of partons, the parton momentum spectra 
in x are mapped almost unaltered onto the hadron distributions. We also conclude, that 
the comparisons in Fig. 13 do not indicate any clear preference for one value of B or L c 
over the other. However, as we will discuss in Sec. 5.3, an indication may be drawn from 
Bose-Einstein correlations among produced hadrons. 

An example of the composition of the final hadron yield in terms of different particle 
species is given in Table 2, where we compare the result of our simulation of e + e~ anni- 
hilation events at Q = 34 GeV with measured particle multiplicities |^]. The remarkable 
agreement is in accord with the presumption that the different particle yields are essen- 
tially determined by the available phase space and the density of hadron states, and not by 
more complex mechanisms. Further experimental constraints by more sensitive measures of 
event shapes such as thrust, sphericity, etc., may be investigated, but it is evident from the 
results shown that our approach yields an overall satisfactory description that withstands 
confrontation with experiment, and encourages us to study more complex reactions in the 
near future. 

5.3 The Bose-Einstein effect 

From Figs. 10 and 11 we can conclude that the distribution of formed clusters clearly 
resembles the picture of preconfinement p7|], which is the tendency of partons produced 
during the cascade evolution to arrange themselves in color-singlet clusters with limited 
extension in both position and momentum space. Since the clusters are the basic units 
within which the final-state hadrons arise, the ensemble of clusters in phase space, as it 
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builds up with time, can be interpreted as a particle emission source with a space-time 
distribution that is determined by the preceding parton evolution. This notion allows us to 



directly relate the dynamics of cluster formation to the well-known Bose-Einstein effect [48], 
which corresponds to an enhancement in the production rate of identical bosons (in our case 
mainly pions) emitted from similar regions in space and time, arising from the imposition 
of Bose symmetry. Enhancements in the mass spectrum of same-sign pion pairs have been 
seen clearly in e + e~-data (for a review see e.g. |44|] ). Let us briefly recall that the Bose 
symmetry imposed on the production amplitude of identical particles from a distribution 
of sources leads to an interference term in the squared amplitude which is only observable 
if the sources are incoherent. Prom the analysis of e + e~-data one finds [Q] that the Bose- 
Einstein effect is reasonably described by a spherically-symmetric space-time distribution 
of sources with Gaussian form 

p(r) = p(0) exp f-^sj > (89) 

where a p is a radius parameter. Such a source leads to an enhancement due to interference 
caused by the identical-particle effect, relative to the rate with no interference, 

b(q) = 1 + X p exp(-a 2 p q 2 ) , (90) 

where q 2 = m 2 n — 4m^, and = \J [p\ + P2) 2 is the invariant mass of the emitted pion 
pair. The degree of incoherence of the source is measured by X p (=1 for complete incoher- 
ence, and =0 for complete coherence), a p measures the source size in fm, or alternatively, 
Op 1 measures the range of enhancement with respect to q in GeV. 

To observe an enhancement in q amongst identical particles, one must compare the 
particle distributions with corresponding spectra in the complete absence of Bose symme- 
try. Thus, in order to get an estimate of the magnitude of enhancement implied by our 
hadronization picture, we proceed as follows. First we evaluate the pion distributions re- 
sulting from a simulation which does not include the Bose-Einstein effect. Then we repeat 
the calculation, but now imposing Bose symmetrization on same-sign pion pairs by assum- 
ing complete incoherence, corresponding to X p = 1. Here we use the method of Sjostrand 



[49 1, which simulates the enhancement due to the identical-particle effect. Finally, we com- 



pute the ratio b(q) of the pion distributions of same-sign pairs with Bose symmetrization 
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to the one without, as a function of the invariant mass q. In Fig. 14a we show the re- 
sulting enhancement bi c {q) for our two previously-used values of L c , confronted with the 



corresponding distribution obtained by the OPAL collaboration [43| at Q = 91 GeV. It 
is remarkable how well this comparison with the experimental data allows us to separate 
the two choices of L c in our calculations. Clearly the value L c = 0.8 fm appears to be 
strongly favored. In fact, the average source size in this case turns out a p = 0.84 fm, which 
is almost identical to the average cluster size determined from Figs. 10 and 11. On the 
other hand, this value is well in the range of the pion source radius determined by OPAL, 
a e p x P = 0.93 ± 0.17 fm with X e p xp = 0.87 ± 0.14. 

We may thus conclude that our presumed identification of L c ~ a p indeed has physical 
relevance that provides a unique relation between the parameter L c and the experimentally- 
observed pion emission source radius. With this important insight, it would be interesting 
to investigate this issue in more detail, because it provides a promising method to extract 
properties of the partons' space-time evolution and cluster formation from the measured 
particle distributions. Since the structure of the perturbative parton cascade development 
is projected locally onto the cluster distribution, which itself maps on the hadron spec- 
tra, the characteristic shape of the Bose enhancement b(q) will depend only on the local 
environment, which may in turn depend on the physical situation (vacuum, as considered 
here, or medium, as, e.g., in deep-inelastic lepton-nucleus scattering or nucleus-nucleus 
collisions). Thus, by comparing, for instance, the Bose-Einstein correlations measured in 
e + e~ — > hadrons to high-energy heavy-ion collisions, one might extract specific features 
of perturbative QCD in a finite-density and -temperature medium, which are absent in 
vacuum |24], 

As an illustrative example of such a comparison, we show in Fig. 14b the ratio 
bo.6f m (Q)/bo.8fm(q) of the curves in Fig. 14a. Although the individual curves in Fig. 
14a are very similar to each other, their ratio is a very sensitive quantity that filters out 
clearly their subtle difference. It is evident that a smaller L c gives rise to a significantly 
stronger enhancement at low masses k ^300 MeV, peaked at about 1.5 times the pion 
mass. Fig. 14b also shows that the results for Q = 34 GeV and Q = 91 GeV are identical, 
even for this sensitive ratio, which implies that the specifics of the Bose-Einstein effect 
and Bose enhancement are independent of the energy, in agreement with what is observed 
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experimentally [44| 



6. SUMMARY AND PERSPECTIVES 

In conclusion, we have presented a novel approach to the dynamics of parton-hadron 
conversion and confinement, based on an effective QCD field theory and a kinetic multi- 
particle description in real time and complete phase space. Our formulation provides an 
extension of the well-understood perturbative QCD parton evolution to account for the 
full space-time history traced from parton cascade development, via cluster formation and 
decay, all the way to the production of final hadrons. The essential points in our approach 
may be summarized as follows. 

(i) We have constructed a scale-dependent Lagrangian that incorporates both parton and 
hadron degrees of freedom. It is manifestly gauge- and Lorentz- invariant, and consistent 
with the scale and chiral symmetry properties of QCD. The introduction of the scale L(r) 
determines locally which are the relevant degrees of freedom around a given space-time 
point r. 

(ii) The formulation recovers QCD with its symmetry properties at short space-time dis- 
tances, and merges into an effective low-energy description of hadronic degrees of freedom 
at large distances. In between the two regimes it interpolates as determined by the scale- 
(L-)changing dynamics, and results in a transformation from partonic to hadronic degrees 
of freedom. 

(iii) The dynamics is described by a set of coupled kinetic equations that derive from the 
field equations of motion, and yield a real-time description in both position and momentum 
space, constrained by the uncertainty principle. 

As a test application, we have considered the prototype reaction e + e~ — ► hadrons where 
the fragmentation of parton jets and their hadronization serves as a generic process that can 
also be imagined as an integral part of more complex reactions. We investigated in detail the 
specifics of the time evolution of parton shower, cluster formation and hadron production 
in phase space, which extends the usual QCD evolution techniques that are limited to 
momentum space and integrated over time. The consistency with experimental data was 
tested, and we found good agreement with measured hadron spectra. A prospective method 
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to extract the characteristics of the space-time development from Bose-Einstein correlations 
of identical hadrons was suggested. Our main results are: 

(i) The details of the parton-hadron conversion are controlled by the quantity L c , the spatial 
separation of neighbouring color charges in their restframe, which defines the scale at which 
non-perturbative confinement forces become substantial. The value of L c is determined by 
the choice of the bag constant B and the condensate value xo- All other macroscopic 
quantities which characterize the evolution, such as pressure, energy density, the time 
scale of the transition, etc., are then determined self-consistently. The found values are in 
agreement with common phenomenology. 

(ii) The time scale obtained for the transition is a remarkable result: the local conversion 
of partons to clusters occurs very rapidly (« 0.1 — 0.2 fm), but the global time scale for 
the transition of the system as a whole is long (« 10 — 30 fm). 

(iii) The QCD features of the perturbative parton evolution are projected unscathed onto 
cluster and hadron distributions, because the conversion is a local, universal mechanism. 
As a consequence the multiplicities of charged hadrons and their momentum spectra are 
predetermined by the preceding parton evolution, which we find in good agreement with 
experiments. 

(iv) Our main result is the sensitivity of the Bose-Einstein correlations among identical 
pions due to Bose symmetry. It allows us to identify the parameter L c with the hadron 
emission source radius measured in experiments, and to fix its value rather precisely to 
L c ~ 0.8 fm. Moreover, the ratio of Bose enhancement of same-sign pions in different 
scenarios or physical situations can provide a very sensitive probe of the environment in 
which the parton system evolves. It can be exploited to study, e.g., modifications to parton 
evolution in finite-density media. 

It is interesting to note that our model appears to be able to correlate sucessfully such 
diverse quantities as the macroscopic bag constant B, the microscopic length scale L c for 
the parton-hadron transition, the associated critical temperature T c , and the magnitude of 
measured Bose-Einstein correlations. 

Finally we comment on future applications. These are manifold, as the advocated 
picture of parton-hadron conversion is universally applicable to any dynamical process 
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where the issue of hadronization arises. Since the strength of our statistical real-time 
description lies in resolving the details of the space-time structure, situations where parton 
cascades undergo interactions with an environment would be the interesting to investigate. 
Let us give three examples: 

a) In deep-inelastic lepton-nucleus scattering the primary quark struck by the photon can 
travel and reinteract before hadronizing, and produce a cascade of secondary partons that 
differs from a parton shower in vacuum. The secondaries are themselves potential candi- 
dates for hard re- interactions, and can lead to a specific A (atomic number) dependence for 
the final-state hadron production. Clearly, here it is essential to keep track of the parton- 
hadron conversion at each point in time and space, because partons that reinteract will not 
be able to hadronize before they approach the free-streaming regime. 



b) In high-energy nucleus-nucleus collisions [51] (y/s/A£ J 200 GeV) the parton density of 
the highly Lorentz-contracted nuclei is very large already in the initial state, and is further 



increased by the materialization and multiplication of partons [24]. Therefore multiple 
scatterings of partons can easily lead to a large number of simultanously-evolving cascades 
that also can interact with each other. In order to resolve such an intertwined structure 
of parton interactions, the space-time dynamics of the system must be taken into account. 
Again, here a microscopic space-time description of parton-hadron conversion is crucial to 
resolve the details of such an intertwined structure of parton interactions and the following 
hadron formation process, depending on the local densities of surrounding parton and 
hadron matter. 

c) The QCD phase transition from a hot, deconfined quark-gluon plasma to excited hadron 
matter as occurred in the early Universe |J7| is of long-standing theoretical interest. Lat- 
tice QCD calculations to date can only investigate the critical behaviour in the vicinity 
of the transition temperature. Moreover, a dynamical evolution of the system deviating 
significantly from thermal equilibrium is not achievable. In view of the future experimen- 
tal programs at RHIC and LHC, it will soon become possible to recreate the QCD phase 
transition in the laboratory |p2| ) and to investigate its dynamics in the real world. It is 
clear that the conversion of a quark-gluon plasma into hadrons is a much more complicated 
process than the hadronization of final-state partons in free space (as in e + e~ annihilation), 
or dilute systems (as in hadron-hadron collisions). In heavy- ion collisions the transition is 
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expected to proceed as a complex evolution of expansion and cooling of the plasma, perhaps 



through a mixed parton-hadron phase, until a purely hadronic phase is reached [53]. There- 
fore, a fully dynamical and microscopic hadronization scheme as proposed in this paper is 
needed to trace the space-time-dependent cooling and expansion process from the parton 
to the hadron phase. 
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TABLE CAPTIONS 



Table 1: Obtained values for phenomenological quantities involved in our description. The 
bag constant B and the condensate value xo = (0| x| 0) are fixed inputs. The values of the 
length scales L c , L x , the global conversion time scale To oc Lq, and the magnitude of the 
surface tension cr c , are then extracted from the numerical simulation. Also listed is the 
associated critical temperature T c , the mass scale of the lightest scalar glueball m x , and 
the estimate for the gluon condensate Go- 

Table 2: e + e~ — * hadrons at Q = 34 GeV: Average multiplicities of partons (n qg ) = 
{n g ) + J2f( n q + n q)-> °f clusters (n c /), and of charged hadrons (n c ^), plus the contribution 



of pions, kaons and protons, in comparison with measured particle multiplicities [47]. 



FIGURE CAPTIONS 

Figure 1: Schematic behaviour of kl(x), eq. (|l~6|) and //(x), eq. ([HI). 

Figure 2: Form of the scale-dependent potential V(L), eq. ([H]), where L x characterizes 
the point of inflection at x x , L c marks the point when the two minima are degenerate at 
Xo and Lq when the potential has a single absolute minimum at xo- The value at x = is 
equal to the vacuum pressure (bag constant) B. 

Figure 3: Diagrammatic representation of a) the two-point Green functions S, D, A, A; b) 
the self energies S, n, 5, H; c) the corresponding Dyson- Schwinger equations (|30|). Notation: 
Fat (thin) lines indicate fully-dressed (bare) propagators, shaded circles and boxes denote 
the full quark and gluon vertex functions, black circles and boxes with attached loops 
represent the local interactions with the collective fields x an d U via the potential V, eq. 

(ID- 

Figure 4: Illustration of the simultanous evolution in space-time r and — 'orthogonal' to 
it — in energy- momentum p of the Wigner functions F a (r,p), according to the transport 
and constraint equations (|||) and (j35|). The self-consistent solution of these equations 
corresponds to summing over all possible quantum paths r, accounting for fluctuations in 
p, under the constraint of the uncertainty principle. 
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Figure 5: Schematics of the different stages of the process e + e~ — > hadrons. The initial 
qq pair with large invariant mass Q initiates a shower which can be followed perturbatively 
until L ~ L x . At this point the conversion into clusters sets in, which is completed around 
L = Lq and followed by the decay of clusters into hadron states. 

Figure 6: Diagrammatic representation of the kinetic equations (|49|). a) The operator 
JC describes free propagation plus the effect of the mean field, b) The integral operators 
A, B, C, . . . include the squared amplitudes for the various interaction processes among the 
different particle species, which change of the particle distributions according to the balance 
of gain (+) and loss (— ) terms. 

Figure 7: a) Space-time evolution of the parton density profile (in arbitrary units) in the 
( r z, Tj_)-plane at different times in the center-of-mass frame of the initial dijet system with 
energy Q = 100 GeV. b) Corresponding development of the cluster density profile as it 
builds up in time due to the conversion of partons. 

Figure 8: a) Time evolution of the kinetic pressures P qg of partons and P x of pre-hadronic 
clusters for gg-initiated jet evolution, for total jet energies Q = 10 GeV (top) and Q = 
100 GeV (bottom). The dashed and full lines correspond to the two parameter choices 
(B 1 / 4 ,xo) =(240,200) MeV and (B 1 / 4 , X o) =(180,100) MeV, resulting in L c = 0.6 fm and 
L c = 0.8 fm, respectively, b) As a), but for a (/(/-initiated jet evolution. 

Figure 9: Total transverse momentum p±(t), eq. (|88|), generated during the time evolution 
of the system in the center-of-mass of the initial dijet system, in correspondence to Fig. 8: 
a) case of (/g-initiated jet evolution (top); b) case of a (/(/-initiated jet evolution (bottom). 

Figure 10: Cluster spectra for L c = 0.6 fm (top), and L c = 0.8 fm (bottom), and total 
jet energies Q = 10 (100) GeV. a) Distribution of the cluster sizes of clusters formed from 
neighboring partons. b) Associated cluster mass spectrum. 

Figure 11: As Fig. 10, but now with the additional constraint of a maximum allowed 
invariant mass per cluster of M cr n = 4 GeV: a) cluster size distribution; b) cluster mass 
spectrum. 
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Figure 12: a) Total gluon and quark multiplcities (n qg ) = {n g ) + J2f( n q J <~ n q)j f° r L c = 0.6 
(0.8) fm as a function of energy Q. b) The corresponding ratios of charged hadrons to 
partons (n c h) / (n qg ) , and of clusters to partons, (n c i) / {n qg ) . 

Figure 13: a) Calculated average charged multiplicity versus total energy Q in e + e~ 
annihilation events, in comparison with experimental data PH| . b) Momentum spectra of 
charged hadrons with respect to the variable ln(l/x), where x = 2E/Q, at Q = 34 GeV 
and Q = 91 GeV, confronted with distributions measured at PEP/PETRA and LEP [46|. 



Figure 14: a) Simulated Bose-Einstein enhancement bi c (q) as a function of the pair 
mass q of same-sign pion pairs for the two values of L c at total energy Q = 91 GeV. The 



data points are from the OPAL experiment [43] at LEP. b) Ratios of the enhancements 



^0.6 /m(<z)/fro.8 fm{q) for total jet energies Q = 34 GeV and Q = 91 GeV. 
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= 230 MeV 
Xo = 200 MeV 


B 1 / 4 = 180 MeV 
Xo = 100 MeV 


L c (fm) 


0.6 


0.8 


L x (fm) 


0.45 


0.65 


t (fm) 


8.5 (21) 


10 (26) 


al /3 (MeV) 


40 


48 


T c (MeV) 


160 


125 


m x (GeV) 


1.05 


1.30 


G (GeV/fm 3 ) 


1.25 


0.50 



Table 1 



Q = 34 GeV 


L c = 0.6 fm 


L c = 0.8 fm 
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